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Abstract 

A simple version of exact finite dimensional reduction for the varia- 
tional setting of mechanical systems is presented. It is worked out 
by means of a thorough global version of the implicit function theorem 
for monotone operators. Moreover, the Hessian of the reduced function 
preserves all the relevant information of the original one, by Schur's 
complement, which spontaneously appears in this context. Finally, the 
results are straightforwardly extended to the case of a Dirichlet prob- 
lem on a bounded domain. 
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1 Introduction 

Conservative classical mechanics can be regarded as a collection of variational 
problems: a mechanical path will be a curve which makes stationary the 
action functional. Such a critical path will be a local minimum for the action 
in the uniform topology in the space of the admissible (w.r.t. the boundary 
conditions) paths if there are no conjugate points along it; this condition is 
known to be related to the positive-definiteness of the second variation -the 
Hessian- of the action functional. A precursor was Jacobi who proposed an 
early theory of the second variation in 19th century. But the very notion of 
conjugate point first appeared in Bliss [3] -see e.g. the historical notes by 
Caratheodory [6] and Eraser |21) . The definitive settlement lays on the well 
known framework of the Morse theory, |26| \25\ l20j . 
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In practice, one is often forced to solve the variational problem numer- 
ically. It turns out that under certain hypotheses it is sufficient to solve a 
reduced problem in a subspace of the space of admissible paths, involving 
a finite number of equations, and recover the very same solution of the full 
problem. Such a technique is sometimes referred to as an exact finite reduc- 
tion method. An early approch of this kind can be recognized in [24j where it 
was employed in the framework of infinite dimensional Morse theory. Later, 
two different reduction schemes have been proposed, both initially used to 
give new versions and proofs for the 'Poincare last geometric theorem': the 
broken geodesic method by Chaperon |1H \T2\ [T3| [T^ [T5] and the Amann- 
Conley-Zehnder reduction [Tl|T6l[T7]- We will focus on the latter approach, 
a concise exposition of which can be also found in |29) . Such a method can 
be employed in many systems, such as fluid [8] and molecular dynamics |28| 
and field theory [7]. 

We propose a simple exact reduction theorem for mechanical Lagrangians 
with fixed endpoints as boundary conditions, giving an estimate of the opti- 
mal dimension of the reduced space (the cutoff N) and comparing our results 
for the cutoff with the ones obtained by means of fixed point techniques 
We then show that the Morse index of the solution can be easily expressed in 
terms of numerically computable quantities; here, as in many other situations 
|10| I18|. the linear algebra trick known as 'Schur complement' spontaneously 
appears. The results are straightforwardly extended -Sect. 5- to the case 
of Dirichlet problem on a bounded il. C M™. 

2 Setting 

We consider a mechanical Lagrangian, 

C{q,q) = ^\e-V{q) (1) 

where the potential energy V € C'^(M",M) satisfies 

sup \V"{q)\=C <oo (2) 

SO that the Hessian of the Lagrangian is bounded: 

sup |Hess[£](g, g)| = max{l,C} = C (3) 

{g,(j)eR2n 

One can now construct the action functional acting on the curves of the 
Sobolev space with fixed endpoints: let 

r = {7 G H\[0,T],R^), 7(0) = go, l{T) = qr} 
L:T3-f{-)^ [ C{j{t),^{t))dt eR 
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Observe that the afHne space F can be parametrized 

lit) = qo + ^^Y^t + cit), cen^H^{[0,T],R") (4) 

where V. is the 'loop space based at 0' (c(0) = c{T) = 0). We define the 
action functional L = L o t on T-L: 

L{c) = L{i{c)), i{c){t) = qo+ ^^^t + c(t) 

It is easy to see that L is twice continuously differentiable in the Frechet 
sense. First we have, for every c,h£T-L 

dL{c)[h] = (VL(c), /i) = ^ (c(t) • h{t) - V\L{c{t)))[h{t)]^ dt (5) 

where for every c £ Ti we can interpret dL{c) € %* as an element VL[c) € % 
by Riesz representation theorem. Next, we find that the 'second variation' 
is the bilinear form in h, k E T-L: 

d^L{c) [h, k] = (^h{t) ■ k{t) - V"{L{c{t))) [h{t), k{t)]^ dt (6) 

which can also be considered as a self-adjoint linear operator on the Hilbert 
space H. We represent every c G H by means of the orthogonal basis 
{fk{'t))k>i,j€{i,...,n}j where for each A; > and j G {!,..., n} the func- 
tion (fj^ : [0, T] —7- has all components 0, except for the j— th which is 
a/2/T sin(fc t/T); to simplify notation, we set (f^ = {fk^---i'Pk)^ ^^^^ ~ 
(c^, . . . , c^) and c^*^) • (/Pfe = X]j=i '^j^ so that we may write 

c(i) = j;c(^').^.(t); ||c||^.:= f\c{ttdt = Y,^,k^\c^'^\'; (7) 
fe>0 •'^ k>0 

notice that then \c(t)\'^ dt = Yjk>Q 

3 Exact finite reduction 

For any given natural number N, using the ordered basis (p'j^, we can decom- 
pose H in two orthogonal subspaces, one of which is finite dimensional (of 
dimension nN): 

n=UN(B Vjv, Un = sp&n{^pi, l<k<N,j e{l,...,n}} 

Adopting a common notation for exact reduction techniques, we define the 
orthogonal projectors PatjQ^v onto Un and V^, respectively. W^r will be 
identified with the subspace of Ti* consisting of functionals vanishing on Vat, 
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and similarly Vj^. When considering the differential of e.g. L on Z/Z^r © VjVj 
we will write dL = duj^L + d\;j^L with duj^L G and d\;^L € V^. By Riesz 
theorem they will be identified with PjvVL and QjvVL respectively. When 
interpreting d^L as a symmetric linear self-operator of we will write 

meaning d^^L = F^d'^LFisf, duj^dy^^L = Q^d'^LF^, and so on. We wih 
also denote by the same symbol operators such as d^^L on H and the one 
that it induces on Un- 

In the following we use [_x\ to denote the largest integer not greater than 
the real number x. 

Proposition 1. For N > and any fixed u G Un, the restricted 

nonlinear operator 

F{u, •) = QjvVL(« + •) : Vjv ^ Vjv 
is strongly monotone. 

Proof. By strongly monotone we mean that there exists some /x^ = /x > 
such that for any vi,V2 G Vat 

(QivVL(u + V2)- QivVL(u + vi), 1^2 - t^i) > ||t;2 - (8) 

We have that 

{QnVL{u + V2) - QnVL{u + -ui), ^2 - vi) = 

= / d^ L{u + vi + X{v2 - vi))[v2 - vi,V2 - vi]dX> 
Jo 

>[ [ (\v2-vi\'^ - sup\V"{q)\\v2-vi\AdtdX, 
Jo Jo \ oeiR" / 



C] \v^^^\\ 



by setting v = V2 — vi, 



k>N 

We wish to prove that there exist a positive integer N and n > such that, 
for every v G Vat we have 



i:(=:i-c)i.«')p>Mii*.=Mi: ^ 

k>N ^ ^ k>N 

this is plainly impossible if for some k > N we have ir'^k^/T^ — C < 0. If 
N =[T VC/tt\ , we see that (tt k/Tf - C > for every integer k> N; let 

/X = mm {{{nk/Tf - C)/(7r fc/r)^} = 1 - CT^i^ {N + l)f > 0; 

A:>iV 

clearly, and fi are as required. □ 
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The cutoff N = [Tv^/vrJ measures, in a sense, the minimum amount of 
information needed to describe the solution of the variational problem, and 
also provides an upper bound to the degeneracy of the number of solutions, 
as we shall see later. Hereafter will be at least as large as [T^/C/ir]. 

From the above calculation it also easily follows that <9y^L is a uniformly 
positive definite self adjoint linear operator on Vn and thus a Hilbert spaces 
isomorphism. 

We are now ready to prove the key statement for the exact finite reduc- 
tion: 

Proposition 2. There exists a continuously Frechet differentiahle map v : 
Un ^ Vat such that, for any u £ Un, 

QN'^Liu + viu)) = 0, or {V L {u + v) , 5v) = V(5v G Vat <^ v = v{u) 

(9) 

Proof. Strong monotonicity of the mapping F{u, •) : Vn — ^ Vat implies that 
this map is a self-homeomorphism of Vat (see e.g. Theorem 11.2 in Chap. 3 
of |19j ; but stronger hypotheses allow us a much simpler proof than the one 
given there, see Appendix). If we rewrite ([9]) as F{u,v) = 0, we immediately 
have that for each fixed u G Un there is a unique solution v{u) G Vn- That 
u I—)- v{u) is continuously differentiahle follows from the implicit function 
theorem: as remarked above d^^L is uniformly positive definite and thus a 
Hilbert spaces isomorphism; and it is immediate to see that this linear map 
is also the partial differential dyj^F{u,v) on the subspace Vn- D 

Remark 1. To obtain an explicit expression for v one can proceed as follows: 
first one has to estimate C < oo, and find the corresponding value of = 
[T-v/C/ttJ as obtained in Proposition[J\ Then by Newton's method one can 
solve recursively for each fixed u 

l)n+l{u) = Vn{u) - [dv,^F{u,Vn{u))Y^ ■ F{u,Vn{u)) 

which converges to the unique fixed point v{u). Observe that the inverse 
d^^F~^ is a bounded operator on Vn\ furthermore, the convergence is faster 
when d\;j^F is far from zero, i.e. when one choses a large A'^ (this also hap- 
pens in [28]). In practice to solve numerically the recursion equation it will 
be necessary to consider a truncation of Vn- One can avoid the inversion of 
dvf^F by suitable approximate methods, e.g. [5} 122). 

Remark 2. Once the function v is known, solving fully the variational prob- 
lem is equivalent to finding all the roots u of the finite dimensional system 
of equations 

PAfVL(ii + {i(n)) = 0, or {V L [u + v{u)) ,5u) = Q ^du gUn (10) 
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By a straightforward argument -see (|16p below- we see that these equations 
are precisely the comeout from stationarizing the reduced action functional 
S: 

^W = o (u) 

where 

Un ='i^^ 3 S{u) = L{u + v{u)) £R (12) 

These equations can be written explicitly in term of the Fourier components 
of U = (Ufc)i<fc<Ar,je{i, ...,„} 

dS vr^A;^ • 

^ = ^ —^ui = [VV{L{u + v{u)))]i, 1< A; <iV,jG {!,..., n} 

(13) 

Finally, to any root u of (fTTj) . there corresponds the admissible mechan- 
ical curve 7 = l{u + v{u)). 



We conclude this section by comparing the minimal size of the cutoff A'^ 
with the known results obtained by other techniques; in particular, in |9] it 
has been shown by a fixed point argument that (in a more general case) it 
is sufficient to take such that 

'l + V2A^)=a<l (14) 



2itN 

It is easy to prove that if verifies the preceding inequality then it is larger 
than [T VC/tt] : in fact, observing that since TV > 1 we have {l + V 2N)/2 > 
1, and recalling that C = max{l, C} 

^A,>ci±|^>c>v^>vc, 

hence 

N >tVC/tt > ItVC/tt\. 

Thus our estimate of the cutoff improves the preceding one, at least 
in the mechanical case. 



4 Morse index 

For a Lagrangian which is positive definite in the velocities, the Morse index, 
being the maximal dimension of a negative subspace of (PL{c), is finite and 
equal to the number of conjugate points counted with their multiplicity [261 
[25l [20] along the critical path 7 = i(c),c = n + v{u), §f ("u) = 0. Recall 
that the conjugate points are the obstruction for such a path to be a local 
minimum, as mentioned in the introduction. 

We want to relate this number to the index of the 'finite Hessian' (PS. 
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Proposition 3. The Hessian of the reduced functional S{u) = L{u + v{u)) 
on Un is 

d^S{u) = dl^L{c) - dv,du,L{c) ■ [dl^m] ■ du,dv,L{c) (15) 
where c = u + v{u) . 
Proof. Differentiating S once, we liave 

dv 

dS{u) = duf,L{u + v{u)) + dvj^L{u + v{u)) ■ —{u) = duj^L{u + v{u)) (16) 

because d\;^L{u + v{u)) = identically by definition of v. Differentiating 
this last identity we get 

= duj,dvj,L{u + v{u)) + dlj^L{u + v{u)) ■ 



which, since L{u + v{u)) is strictly positive, gives 



dv 
du 



Differentiating (11611 once more, we get 



{u) = -[dl^L{u + v{u))\ ^ -du^dv^Liu + viu)) (17) 



d^S{u) = dl^L(u + i(u)) + dv^du^Liu + v(u)) ■ —(u) 



that together with (jl7p gives the thesis. 

The Hessian, on T-L = IAn ® Vat, can be written in the block form 



□ 



9u^L dv^du^L 



where (du^Ovf^L)* = dvf,du^L. 

Proposition 4. The Morse index iL{u + v{u)) of a solution equals the index 
of S at u: 

iL{u + v{u)) = is{u) (18) 
The same is true for the nullity. 

Proof. A standard technique (see e.g.|10 | [T8 ] ) to block-diagonalize the bilin- 
ear form (f^L is by means of the block matrix 



T 



-1 



which clearly defines a bounded linear operator on H = Un © V^v, with 
bounded inverse 



The same holds for the adjoint T* . We can thus use them to perform a 
change of basis on so that dpL takes the block-diagonal form 

-1 



T*(fLT 



( rP- 



V 



It follows immediately that negative and null spaces of (PL{u + v{u)) and 
d'^S{u) are isomorphic via T. □ 

Remarks. Contained in the previous argument is the a priori upper bound 
for the index: 

iL < dimWjv = nN (19) 



5 Dirichlet problem 

Generalizing the preceding one dimensional case, wc consider an open bounded 
Q C M*" with regular boundary di^; in this case the Hilbert space we need 
is n = H^in,R), the closure of in H^{n). The variational problem 
VL(J)) = 0, ^ G for the functional 



m = j^i^^ wm\' - v{<t>{x)) j dx (20) 

can be regarded as a weak formulation of the Dirichlet problem 

A</,= -y'(0), <t>\9n = 

In this case the basis vectors will be the eigenfunctions of the Laplacian; as 
customary we arrange the distinct eigenvalues of the opposite of the Lapla- 
cian in increasing order, with each Afe repeated according to multiplicity: 
< Ai < A2 < A3 . . . , and A^ — t- +00 as A; — t- 00; an orthogonal basis will 
be {ipk '■ k > 1}, where —Aipk = Xkfk and ^Pklea = 0. Of course all this 
depends on Q. For every integer > 1 we define 

Un = span{<fk eH = H^{n,R) : k < N} = FnH 

and Vat = QnH is the orthogonal complement of Un- Assuming that 

sup|y"((/))| =C<+oo (21) 
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we can again write (recalling that |Vf 




{QnVL{u + V2) - QnVL{u + vi),V2 - vi) > 



> [ {\V{v2-vi)\^ -C\v2-vi\^)dx= [ {\Vv\^-C\v\^) 




= 5: (A. - C) l.^'^) 



k>N 



Pick the smallest integer N such that Xn+i > C and set 



= min{(Afc - C)/Xk : k > N} {= I - C/Xn+i) ; 



we have that fj, > and as before 



k>N k>N 



giving the required strong monotonicity. In particular, one gets that the 
Morse index ^^(c) of any solution c = u + v{u) equals its reduced index 
isiu), which is always bounded from above by the dimension of Un. 

Remark 4. Recall that for m = 1 and vector valued functions in we found 
that the minimal dimension of to achieve the reduction was nN, with 
N = [T \fC I ttJ . This dimension is exactly the number of the eigenvalues of 
—A not larger than C; when m > 1 it grows as with C — )• 00 as stated 



where is the Euclidean ball in M™. Let us remark that the Morse index, 
even in the case m > 1, can be related to the multiplicity of the solutions 
via the min-max theory (see e.g. [271 [3]). Furthermore the solutions and 
their geometrical features can also be studied. In fact, for a given solution 
define the nodal set Z C as the closure of {x € : <^{x) =0}. It can be 
shown [30j that the Morse index provides an upper bound to the norm 
of the solution, to the Hausdorff measure of Z and to the vanishing order of 
X G Z; in our case this provides an a priori estimate of such properties of the 
solutions, since the upper bound on i^. depends only on C = sup^gjg |1^"(0)| 
and on the domain fi. 

Appendix 

We present here a simple, self contained argument, for the proof of Proposi- 
tion 2. 



by Weyl's law [23] 



dim(Z^Ar) 



(27r)™ 



(22) 
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Proposition 5. Let Y be a real Hilbert space, X a Banach space, D an 
open subset of X. Assume that F : D x Y ^ Y is a map, and that for 
every x € D the x— section of F, i.e. the map F{x,-) : Y ^ Y given by 
y i-> F{x, y) is strongly monotone, that is, there exists fj.^ > such that 

{F{x,u) — F{x,v),u — v) > fix \\u — vW"^ for every u,v €Y. 

Then the zero-set Z{F) = {(x,y) & D x Y : F(x,y) = 0} is the graph of a 
map if : D ^ Y. 

Proof. The point x £ D being given, denote for simplicity by / the j;— section 
of F, that is, f{u) = F{x,u) for every u £ Y, and let /i = It is 

immediate to see that for all pairs u,v E y we have ||/(ii)— > jj, \\u—v\\. 
This clearly shows that / is injective with a Lipschitz continuous inverse 
g : f(Y) — )• Y; moreover / maps closed subsets of Y into closed subsets of 
Y: in fact if Vk = f{uk) converges to v G Y, then Uk = g{vk) is a Cauchy 
sequence in Y; if u is its limit, then v = f{u) by continuity of /. In particular, 
fiY) is closed in Y . Next observe that for every c GY the differential df{c) 
of / at c is also strongly monotone (with the same constant /x). In fact, given 
/i G y we have, for every t G M: 

{f{c + th) - f{c),th) > fx\\thf; 

assuming t non zero: 
and when t — )■ we get 

{df{c)[h],h)>fi\\h\\'' 

as desired. 

By what we have just proved the range (i/(c)[y] is then a closed subspace 
of Y] and if v is orthogonal to it we get 

Q={df{c)[vlv)>fi\\vf, 

so that V = Q. Thus df{c) is onto, and hence an isomorphism of Y onto 
itself. The local inversion theorem is then applicable to /, proving that it 
is a local diffeomorphism at every point of Y] then f{Y) is open, and being 
also closed it coincides with Y . We have proved that for every x £ D the 
a;— section F{x, •) is a self-diffeomorphism of Y. The conclusion is now 
immediate: for every x € D let (p{x) be the unique zero of y i— F{x,y); the 
implicit function theorem (see, e.g. ^) proves that ip is of class C^. □ 
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Remarks. li X = Un, Y = Vn and L = u, in Section 2 we have the 
following setting: a scalar valued function u : X xY ^ M., twice continuously 
differentiable in the Frechet sense, such that the second differential with 
respect to Y is strongly positive, i. e. there is > such that 

^u{x,y)[h,h]>ii\\hf, 

and moreover lim||j,||^oQ u{x,y) = +00 (also implied from this positivity). It 
is well-known that these hypotheses imply that the x— section of u attains its 
absolute minimum in a unique point <p{x) G Y, for every x £ X. The proof of 
this genuinely variational feature is based on the weak compactness of closed 
bounded convex subsets in a Hilbert space; this point is clearly the only 
singular point of the differential 9yu(x,y) G 1"*, that is, the unique point 
ip{x) G Y such that dYu{x, ip{x)) = 0. This implies, fully inside a variational 
scenario, the existence and uniqueness of ip. However, even though this 
argument is undoubtedly elegant, smoothness of ip does not follow, and this 
technique does not seem to be simpler than the one proposed in the previous 
Sections, which involves only elementary analysis, and no weak compactness 
arguments. 

Acknowledgements. We thanks Giulia Antinori, Marco Degiovanni and 
Luigi Salce for fruitful conversations and suggestions. 
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